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Abstract 

It remains to be determined experimentally if massive neutrinos are Majorana or Dirac particles. 
In this connection, it has been recently suggested that the detection of cosmic neutrino background of 
left-handed neutrinos and right-handed antineutrinos in future experiments of neutrino capture 
on beta-decaying nuclei (e.g., + —?► ^He-|-e“ for the PTOLEMY experiment) is likely to distinguish 

between Majorana and Dirac neutrinos, since the capture rate is twice larger in the former case. In this 
paper, we investigate the possible impact of right-handed neutrinos on the capture rate, assuming that 
massive neutrinos are Dirac particles and both right-handed neutrinos and left-handed antineutrinos 
can be efficiently produced in the early Universe. It turns out that the capture rate can be enhanced 
at most by 28% due to the presence of relic and 17^ with a total number density of 95 cm“^, which 
should be compared to the number density 336 cm“^ of cosmic neutrino background. The enhancement 
has actually been limited by the latest cosmological and astrophysical bounds on the effective number 
of neutrino generations = 3.14to;43 at the 95% confidence level. For illustration, two possible 
scenarios have been proposed for thermal production of right-handed neutrinos in the early Universe. 
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1 Introduction 


Although a number of elegant neutrino oscillation experiments in the past few decades have well established 
that neutrinos are massive particles, it is still unclear whether massive neutrinos are of Dirac or Majorana 
nature mm- Thus far, tremendous efforts have been placed on the experimental searches for neutrinoless 
double-beta (Oz//3/3) decays, which take place only if lepton number violation exists and massive neutrinos 
are Majorana particles Oil El [ 6 ]. The experimental discovery of duj3j3 decays will provide us with a robust 
evidence for Majorana neutrinos. However, in case that Oi^/3/3 decays are not detected in all the future 
Oz//3/3 experiments, it is still possible that neutrinos are Majorana particles, if neutrino mass ordering is 
normal (i.e., < m 2 < m 3 ) and an intricate cancellation occurs in the effective neutrino mass relevant 

for Oi^/3/3 decays (see, e.g., Ref. 0)- In this case, another independent approach should be utilized to 
probe the Dirac or Majorana nature of massive neutrinos. 

More than fifty years ago, Weinberg pointed out [ 8 ] that the cosmic neutrino background (Ci/B) 
predicted by the standard Big Bang theory of cosmology can be detected via neutrino capture on beta- 
decaying nuclei, e.g., —)■ ^He -|- e~. This possibility has been extensively studied in many recent 

works [aiiniiiiiiiiiiiaiii!. In particular, for the future experiment PTOLEMY m with 100 grams of 
tritium, the capture rate r(z/g -|- —)• ^He -|- e~) has been found to be [H] 

Tm ~ 8 yr“^ (Majorana) ; T^ Ri 4 yr“^ (Dirac) . (1) 

These results have profound implications for cosmology and elementary particle physics. First, a successful 
detection of Ci^B is very important to further verify the standard theory of cosmology [m dzi US], and 
serves as a unique way to probe our Universe back to the time when it was just one second old. We 
already have an excellent example that the precise measurements of cosmic microwave background (CMB) 
have given valuable information on the Universe at the age of 3.8 x 10^ years, and greatly improved our 
knowledge on the cosmology. Second, the relation Tj^ = 2 r£) between the capture rates in Eq. ([T]) offers 
a novel way to distinguish between Dirac and Majorana neutrinos. In this paper, we concentrate on the 
second point and take it more seriously. 

It is worthwhile to emphasize that the contributions from right-handed components of massive Dirac 
neutrinos are completely neglected in the calculations leading to Eq. ([T]). See, Ref. [H], for more details. 
An immediate question is how the right-handed Dirac neutrinos are produced in our Universe, in the 
standard theories of particle physics and cosmology, and whether their abundance can be safely neglected. 
The second question is how the right-handed Dirac neutrinos affect the detection of CzzB, i.e., the capture 
rate in Eq. ([1]), if they are copiously generated in the early Universe and survive today as a cosmic 
background. In order to answer these two questions, we assume that massive neutrinos are Dirac particles, 
and investigate carefully their production and evolution in the early Universe, both within and beyond 
the standard model of particle physics (SM). 

The remaining part of the present paper is organized as follows. In Sec. 2, the thermal production 
of right-handed neutrinos in the minimal extension of the SM with massive Dirac neutrinos is reviewed. 
The production rate turns out to be extremely small and can be neglected. Then, we investigate the 
cosmological constraint on the relic density of right-handed neutrinos in Sec. 3, assuming that they can 
be thermalized in the early Universe in the scenarios beyond the SM. Subsequently, in Sec. 4, two possible 
scenarios have been presented to show that they can indeed be thermally produced if the primordial 
magnetic fields or secret interactions among right-handed neutrinos exist. Sec. 5 is devoted to the impact 
of relic right-handed neutrinos on the detection of CzzB. Finally, we summarize our main results in Sec. 6 . 
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2 The Extended SM 


We first briefly review the minimal extension of the SM with three right-handed neutrinos (for a = 
which are singlets under the SU(2 )l x U(1)y gauge symmetry. The relevant Lagrangian reads 




^aL 0^i')al3 + ^-C. 


(2) 


where £gj^ stands for the SM Lagrangian, ^ ~ denote respectively lepton 

and Higgs doublets, for a,/3 = e,fj,,T are Dirac neutrino Yukawa couplings. After the Higgs field 

acquires its vacuum expectation value (H) = v ^ 174 GeV and the gauge symmetry is spontaneously 
broken down, one obtains the Dirac neutrino mass matrix M-^ = Y^^v, which can be diagonalized via 
a bi-unitary transformation = diaglm^,m 2 ,m 3 }. In the mass basis, and constitute a 

massive Dirac spinor i/j = In the following, we refer to this minimal extension of the SM with 

massive Dirac neutrinos as the extended SM. As indicated by the latest Planck results on CMB m, the 
sum of neutrino masses S = m^ +m 2 + RT -3 is strictly constrained, i.e., S < 0.23 eV at the 95% confidence 
level (C.L.). Therefore, for Dirac neutrinos, we are left with two serious problems. First, a global U(l) 
symmetry corresponding to the lepton number conservation has to be imposed on the generic Lagrangian 
in Eq. ([2|) in order to forbid a Majorana mass term which otherwise is allowed by the SM gauge 

symmetry. Second, the neutrino Yukawa coupling constants = m^/v ~ 10“^^ in the mass basis are 
smaller by twelve orders of magnitude than the top-quark Yukawa coupling ~ 0{1). This exaggerates 
the fermion mass hierarchy problem of the SM. In a realistic model of Dirac neutrinos, these two problems 
should be properly addressed. However, we temporarily put them aside and focus on the cosmological 
implications in the following discussions. 

After specifying the theoretical framework, we are now in a position to consider the production of 

in the early Universe. In fact, this task has already been accomplished in Ref. [20]. However, it is 
instructive to revisit this problem in view of recent progress in neutrino physics (e.g., the establishment of 
massive neutrinos) and the discovery of Higgs boson (e.g., the observation of Higgs-fermion interactions). 
Following the notations in Ref. m, we can calculate the number density of a particle species a, which 
is involved in the interaction X o o -|- Y, via the Boltzmann equation 


dn 
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where (for i = a,X, Y) are the number densities in thermal equilibrium, and H is the Hubble expansion 
rate. Those two terms in the parentheses on the right-hand side of Eq. Q stand for the absorption and 
production rates, respectively. More explicitly, the collision term is given by 


7 (Y^a + Y) 


d^PY d^p^ d^py 
(27r)^2p^ (27r)^2p^ (27r)^2py 


{2TTfS‘^{px - Pa 


Py)e-Px/r|7W(Y^o + Y)p , 


( 4 ) 


where (for i = a,X,Y) are the four-momenta, T is the temperature, and |A4(Y —>• a -|- Y)P should be 
summed but not averaged over the internal degrees of freedom of the initial and final states. In general, 
X and Y can also be a set of multiple-particle states. 

In the early Universe, when the temperature is extremely high T 3> Tew ~ 200 GeV being 

the critical temperature for electroweak phase transition, the SM gauge symmetry is restored and all the 
SM particles are massless, except for the Higgs boson. The right-handed neutrinos only experience 
the Yukawa interactions, as given in Eq. ([2|), so the production and absorption of r'jE ^-re governed by 
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the tiny Yukawa couplings Note that we are working in the mass basis in the sense that the Yukawa 
coupling matrix is diagonal. In this case, the dominant processes for production should be Higgs 
boson decays H —)■ 17 ^ + top-quark scattering + ^ H ^ z7^ + and gauge boson scattering 

V + V ^ H ^ T7^ + Taking only the decays and inverse decays into account, we can immediately 
figure out the decay rate Tjj = y‘fM^/{32TT), and thus the corresponding collision term 

]\/T '~p‘2 

^^iK) = — K^{Mh/T) = 7d , (5) 

where is the n-th order modified Bessel function. The contributions from the scattering processes are 
on the same order as that from decays. Because of finite-temperature effects [2T], the gauge interactions 
result in thermal lepton masses MJ{T) = (3y^ -|- y'^)T^/32, implying a slight reduction of the decay rate 
compared to the result at zero temperature. Moreover, both gauge and top-quark Yukawa interactions 
give rise to a thermal Higgs mass Mfj{T) = (8M^ + M'^ + 2m‘f + — T|^)/(8u^), where the gauge 

boson, top-quark and Higgs boson masses are evaluated at T = 0. Given the latest values = 173 GeV 
and Mfj = 125 GeV, we can obtain M^{T)/T « 3/4 in the limit of T S> As a consequence, the 

Higgs mass Mjj in Eq. ([5|) should be replaced by the thermal one M^{T) ~ 3T/4. Both Higgs bosons and 
left-handed neutrinos are well in thermal equilibrium due to the efficient gauge interactions, so Eq. ([3|) 
can be simplified to 
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where the in-equilibrium density = T'^ for the Boltzmann distribution with a vanishing chemical 
potential. Eor comparison, the number density of photons in thermal equilibrium is given by = 2n^^j^ = 
2T^/7r^. Now it is evident from Eq. (|6]) that if the creation rate Tj^ ^ = 7 d/^ 4 much smaller than the 
Hubble expansion rate H = 1.66y^T^/Mp[, where y* is the number of relativistic degrees of freedom and 
Y/pj = 1.2 X 10^® GeV is the Planck mass scale, the production of will be inefficient. With the help of 
Eq. dH), it is straightforward to derive 


R = T^ /H 


10 


-3 y? ^Pl 


(7) 


T ’ 

which is valid for T ^ Tew Vi ~ 10“^^ and = 106.75 in the SM at T = lO’^ GeV, one arrives 

at ii ~ 10“^®, which is larger by several orders of magnitude than the result in Ref. [20], where only 
neutral-current interactions are included. For T = 10^ GeV, we have an even larger ratio R ~ 10“^^. 
However, such a small production rate indicates that will never be populated in the early Universe, 
which should hold as well at any temperature above T^^. 

Below the electroweak phase transition, i.e., T <C 7 ew> the Higgs bosons, top quarks and weak gauge 
bosons have already decayed away, and the dominant production channel of is through the conversion 
from ‘iue to the presence of Dirac masses [2U] . Given the production rate of being T^^ ~ G|T^, 

for T = 10 GeV; 


G|T^m^, and thus R 


10 


-10 


where G-p is the Fermi coupling constant, we obtain T 
or i? ~ 10“^^ for T = 1 MeV. Therefore, it becomes clear that in the extended SM can never be 
abundantly produced in our Universe. 


3 Cosmological Constraints 

However, in some new physics scenarios beyond the SM, right-handed Dirac neutrinos can be copiously 
generated in the early Universe. Before going to any details of new physics models, we simply assume 
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that can be thermally produced and then decouple from the plasma of SM particles at a freeze-out 
temperature which needs not to be specified at this moment. With such a setup, our discussions will 
be applicable to more general cases. In addition to the Ci^B, will be a new kind of cosmic background, 
which is restrictively constrained by the cosmological observations, such as Big Bang Nucleosynthesis 
(BBN) and CMB. 

At the freeze-out temperature, where the decoupling is taken to be instantaneous for simplicity, the 
number densities of and are equal, i.e., n^^(T^) = n^^{T^). In the later evolution of the Universe, 
t'L remains in good contact with the other SM particles via gauge interactions, and its number density 
is determined by the Fermi-Dirac distribution with a temperature T. However, gets diluted by the 
expansion, and its number density at temperature T is given by 


nu^{T) 



. a{T) _ 


( 8 ) 


where a(T) is the scale factor as a function of the temperature T. On the other hand, assuming adiabatic 
expansion of the Universe, we find that the conservation of entropy in the thermal bath leads to 


9, 


,{T)T^ 


9*siT^)iT,^r 


a(r) 


-I 3 


(9) 


where stands for the number of effective degrees of freedom. Note that we have = g^,{T) before 

the neutrino decoupling. Since the left-handed neutrino is relativistic and in thermal equilibrium before 
it decouples when T > 1 MeV, its number density scales as oc T^. Therefore, at the temperature 
when begins to freeze out, the ratio of the number densities of and r'L reads 

g^^{Tiy 


where the relations and have been used. The ratio of 

right-handed and left-handed neutrino number densities is fixed by g^,g at the freeze-out temperature of 
since we know « 1 MeV and 5*^(7}^) = 10.75. If we have = 106.75, and thus the 

right-to-left ratio of neutrino number densities is ~ 0.1. Since T < both and are freely 

streaming, the above ratio is then unchanged at the present time. As is well known, the Ci^B consists 
of left-handed neutrinos and right-handed antineutrinos whose number density is 336 cm“^. Now 
we also have right-handed neutrinos ^rid left-handed antineutrinos as a cosmic background with a 
total number density 34 cm“^. Even larger number densities of and V-^ are also possible for a lower 
freeze-out temperature, as we shall show later. 

The presence of extra radiation, such as under consideration, could substantially modify the pre¬ 
dictions for primordial abundances of light nuclear elements by the standard BBN theory and the power 
spectrum of CMB. As usual, the contribution of extra relativistic particles x to the total energy density 
of radiation at the CMB temperature is parameterized as 


Pr P'y T Px 


1 + Apff - 



( 11 ) 


where = 3.046 is expected if only three generations of neutrinos in the SM exist. That the value is 
not exactly three can be ascribed to the non-thermal distortion of neutrino energy spectra, as neutrinos 
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Figure 1: The extra effective number of neutrino generations — 3.046 during the era of Big 

Bang Nucleosynthesis has been calculated by choosing different values of the freeze-out temperature of 
right-handed neutrinos T^. The shaded area represents current cosmological and astrophysical constraints 
ANgfj = 0.101 q 43 at the 95% C.L., where the Planck TT -|- lowP -|- BAO data sets [H] are combined with 
the helium abundance measurements [24]. The number of relativistic degrees of freedom has been 

taken from Refs. [25] [26]. 


are still slightly coupled to the thermal bath when electrons and positrons annihilate into photons 
Now the contribution from should also be taken into account. In a similar way, one can derive the 
relationship between the energy densities of and at the CMB temperature Tcmb ~ follows 

Pu„ (^cmb) 


9*s{Tcub)V^" fllY' 


Pi'A'^CMb) I 9'*s(^fo) -I \ ^ J 


eff 


( 12 ) 


where = 3 because of three generations of neutrinos, and = 3.046 originating from the left-handed 
neutrinos The extra due to the presence of right-handed neutrinos can then be defined as 
ANgfj = Ngfj — Agg, and it is found to be 

in 4/3 

i) 


AAgff = 


9*s{T^ 


CMB^ 


11X4/3 

A" 
4 / 


9.siTA) \ 

Assuming that the left-handed neutrinos decouple instantaneously, the above AAg^ reduces to 

'9.siTty"^" 


(13) 


AAge = 


9*s{TA)\ 


A" 


(14) 


in agreement with the results given in Refs. mm [29] 0 According 


to the latest results from Planck 


^To account for the non-instantaneous decoupling effects of we adopt the value of g*s(TcMB) ~ 3.931 in the following 
calculations, while its instantaneous decoupling limit should be 3.909. 
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( 15 ) 


Collaboration, the effective number of neutrino generations is determined to be m 

TVgff = S.Utoit , He + Planck TT + lowP + BAG , 

at the 95% C.L. Any additional relativistic species present in the cosmic background will be stringently 
constrained by Eq. (I15p . In Fig. [1] we have calculated the extra number of neutrino species by 

varying the freeze-out temperature T^. The cosmological bound is represented by the shaded area, while 
electroweak and QCD phase transitions are indicated by the dashed lines. In our calculations, the values 
of g^siT) have been taken from Refs. |25l I26j . Some comments are in order. First, the ranges of 
below 200 MeV are excluded by cosmological observations. This indicates the importance of QCD phase 
transition in diluting primordial relativistic particles. Second, for ~ 200 MeV, the upper bound on 
AVgg in Eq. (fTC|) can be saturated, namely, AV^g = 0.53. In this situation, it is straightforward to get 
~ 0.28. Or equivalently, the number density of relic and Eg is 95 cm“^. As we show later, 
these results affect significantly the detection of CzzB. 


4 New Physics Scenarios 


Now we propose two possible scenarios to realize a thermal production of in the early Universe, and 
the freeze-out temperature should be above Tgco ~ 200 MeV to evade cosmological bounds. 

Primordial Magnetic Fields .—An important intrinsic property of massive Dirac neutrinos is that they 
can have nonzero magnetic dipole moments |30t I31j. If the SM is extended with massive Dirac neutrinos, 
one can obtain [32l |33] 


l^v 


3 X 10 


-20 


m. 


0.1 eV 


Mb 


(16) 


where /ig = e/2mg is the Bohr magneton. It is very likely that primordial magnetic fields exist in the 
early Universe, e.g., resulted from the electroweak phase transition [3H [35l |36], so the magnetic dipole 
interaction leads to an efficient production of through the spin-flipping process r'-g —)• [37]. On 

the other hand, those primordial magnetic fields could also survive until today and serve as seed fields to 
explain the observed galactic magnetic fields around ~ 10“® G [38]. Even though the magnetic fields 
can be generated during cosmological phase transitions, it remains unclear how the random magnetic 
field fluctuations are transformed into macroscopic-scale magnetic fields |38| ■ For simplicity, we follow the 
phenomenological approach in Ref. m and assume random magnetic fields with a scaling behavior as 

1 2 


B{t, L) — Bq 


_a{t)_ 


(17) 


where the term involving the scale factor a{t) indicates the conservation of magnetic flux B ~ during 
the expansion of the Universe. In addition, the index p accounts for how the field strength depends on 
the physical spatial scale L. The initial domain size Lq and field strength Bq are determined by the 
production mechanism operated at the scale factor Oq. In the radiation-dominated epoch, the relation 
a{t) oc oc T~^ holds, so one can convert the dependence of B on the time into that on the temperature 
T. For the small-scale random magnetic fields S> Eg, the transition probability is [37| 

Fl^r = \pil B^LQH-^L^ , ( 18 ) 

where = T*^* « 30^1 is the total weak interaction rate and the inverse of Hubble expansion rate 
comes in as the largest time scale. For a suitable field strength Bq and domain size Lq, there is no 
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doubt that can be brought into thermal equilibrium with In order to ensure the decoupling of 
at latest around ^QCD ~ 200 MeV, we require the transition rate to be smaller than H at ^QCD> namely, 

/T \V2 

^h) ^ 6.7 X 10 G > (19) 

where ^h(^qcd) = -^~^(^qcd) ~ 3.5x 10^ cm and L^(Tq(-,j^) Ri 1.6x 10“^ cm. Given the prediction of 
in Eq. (llOp . we can translate this constraint into an upper bound on the magnetic field Bq generated from 
the electroweak phase transition at if the domain size fulfills Lg > L™“(r) = 10“^ cm (MeV/T). 

Although the primordial magnetic fields in this case dissipate away at the time of BBN, it will not affect 
the production of that is already completed at T > Taking p = 1/2 for example, we can get 


Bo < G 


'3 X 10-2Vb' 




( 20 ) 


where Lq = T™™(rQQj-,) = 5 x 10“^ cm is input. It is now evident that for Bq « lO^^^ G and Lg > L™™, 
as predicted by a specific model in Ref. [34], can be in thermal equilibrium at earlier times and then 
decouple from the thermal bath just before the QCD phase transition. 

Secret Interactions .—It is also possible to populate via an exotic interaction, which is introduced 
exclusively for right-handed neutrinos and thus named as “secret interaction”. For our purpose, it is 
enough to follow a phenomenological approach and simply add two interaction terms into the SM 


C D - g^xi^xV^ , (21) 

where denotes a light vector boson V (e.g., niy ~ 1 MeV), while y a heavy Dirac fermion (e.g., 

~ 2 TeV) as a candidate for cold dark matter. Here g^, and g^ stand for the coupling constants of V 
with and y, respectively. This kind of secret neutrino interaction was first considered for the SM left- 
handed neutrinos [391 SOI E], and later extended to sterile neutrinos of eV-scale masses |42l 1431 HU 1451146j . 
In these works, the secret interaction also applies to a dark matter particle, offering an intriguing solution 
to the small-scale problems of cold dark matter in the cosmological structure formation [39]. 

In this model, we further postulate that the dark matter sector y is also coupled to the SM via some 
high-energy dynamics, such that y can be thermally produced. However, the decoupling between these 
two sectors also takes place quite early, e.g., at ^ 10 TeV. Then, the new interactions in Eq. (|2T|) will 
bring both V and into thermal equilibrium with y, if the coupling constants are not extremely small. 
The evolution of the whole system can be summarized as follows. 

• The dark matter annihilation yy —VV and yy —)• (negligible for <C g^) will be frozen out 

at leading to a correct relic dark matter density for = 2 TeV and g^. ^ 0.8, as 

demonstrated in Ref. |39j . The Vr-x elastic scattering could keep them in kinetic equilibrium until 
T ~ keV. This salient feature, together with the self-interaction of y mediated by V, could help 
solve all the small-scale problems of structure formation. 


• Since Vr is only in contact with the dark sector (i.e., V and y), which deconples from the SM sector 
at Tg, it just cools down as where represents the temperature of 

the SM sector. At the time of BBN, the extra number of neutrinos can be estimated as [42] 


AA^eff = 


Pu^ + Pv 

[PuJ'i) 


/ 10.75 \ 
V106.75 ) 
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3-F - X - 

2 7 


4/3 


0.22 


( 22 ) 







which is compatible with the BBN constraint. Later on, a reheating of and arises from 
V —>■ when V becomes non-relativistic at T < my. The ratio of the temperatures after and 
before reheating is = (11/7)^^^, whereas the photon temperature is increased by electron- 

positron annihilation, so the effective number of neutrinos in the epoch of CMB turns out to be 




Pu 


(Pu^ / 3 ) 


3 X 


9^s v^CMb) 

. 9*s{Td) 



X 



0.26 , 


(23) 


where 5 * 5 ( 705 ^^ 3 ) = 3.91 is used. Therefore, this model survives both BBN and CMB bounds. The 
number density of 1^3 and Vy at present is 34 x 11/7 ~ 53 cm“^, as expected for a high freeze-out 
temperature and a late-time reheating from V decays. 


It is worth mentioning that an interesting scenario of non-thermal production of 1^3 has been recently 
presented in Ref. m, where the coupling between inflaton and z /3 is introduced to generate a degenerate 
r '3 gas via inflaton decays. After its production, 1/3 evolves separately and is diluted as the Universe cools 
down. Although the energy density of 1^3 is constrained by the BBN and CMB observations, as in our 
two scenarios, the relic number density of 1^3 and 17y can be as high as one half of the photon density, 
i.e., around 220 cm“^. The main reason is that 1^3 is conhned in the low-energy part of the non-thermal 
distribution function, compared to the thermal spectrum. 


5 Impact on the Detection of Ci^B 

Finally, we examine the impact of relic right-handed Dirac neutrinos on the detection of Cz^B in a future 
experiment, such as PTOLEMY m- Experimentally, one studies the spectrum of the emitted electrons 
(e.g., from i/g -|- —>• ^He -|- e“), and looks for events that have kinetic energies above the beta-decay 
endpoint. However, because of tiny neutrino masses and small momenta carried by relic neutrinos, the 
expected signals are very close to the endpoint. Eor instance, if three active neutrinos have a degenerate 
mass ~ ttiq, we then expect the signal to show up at the position that is 2mQ beyond the endpoint. 
An energy resolution comparable to the neutrino mass tuq is then required in order to clearly select the 
signal from the dominant beta-decay background. In the recent proposal PTOLEMY, new techniques are 
suggested to achieve such a high energy resolution. The beta-decaying nuclei will be deposited onto 
some surface substrate so as to reduce the nucleus recoil. An energy resolution A ~ 0.15 eV, defined as 
the full width at half maximum (EWHM) of the Gaussian distribution, is claimed to be achievable, so 
the detection of Ci^B is promising, in particular for nearly-degenerate neutrino masses and relatively large 
values of absolute neutrino masses (e.g., ttiq > A). 

Beside the location of the signal, one can also measure its height, which is related to the capture rate. 
Adopting the calculation from m, we find the capture rate of relic neutrinos 

3 

Tci/B = , (24) 

s,=±l/2 i=l 

where stands for the two helical states, and j indicates three different mass eigenstates. In addition, 
n-{Sj^) represents the number density of incoming relic neutrino in the helical state while Nj, = 
Mrp/m3}i is the number of target nuclei. Approximately, we have (yj{s^)v^ ~ where (for 

j = 1, 2,3) are the matrix elements in the first row of the lepton mixing matrix, and a ~ 3.834 x 10“^^ cm^ 
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for the neutrino capture on tritium. More importantly, the spin-dependent factor A{sj^) is given by |14] 



1 — Vuj, = -1-1/2 right helical 

1 + = —1/2 left helical , 


(25) 


where the neutrino velocity Vi,^ = I Pj, I /Ey approximates to 0 and 1 in the non-relativistic and relativistic 

J i 3 

limits, respectively. According to the standard cosmology, the present temperature of Cz^B is around 
Ty = 0.168 meV. Therefore, for the degenerate mass region that can be probed by PTOLEMY, relic 
neutrinos are highly non-relativistic today. In this case, we have A(±l/2) = 1 for both left and right 
helical states. As shown in Eq. if both 13-^ and 13^ are absent in the cosmic background, we have the 
capture rate Tj-, ~ 4 yr“^ for 100 grams of tritium. This rate will be modified in the following cases. 

1. If both 13^ and are thermally produced and decouple from the thermal bath at a high temperature 
above the electroweak phase transition. The number density of right- and left-helical neutrino states 
will be increased by 10%, implying a capture rate T^ ~ 4.4 yr“^. However, if they freeze out just 
before the QCD phase transition, the capture rate can be enhanced by 28%, namely, Tp ~ 5.1 yr“^. 
The latter scenario saturates the upper bound on the extra effective number of neutrinos in the BBN 
and CMB eras. 

2. Unlike the thermal production, the non-thermal and degenerate gas of i3^ and considered in 
Ref. [47] could change the capture rate by 64%, namely, T^ ~ 6.6 yr“^. This modification will 
diminish the chance to distinguish between Dirac and Majorana neutrinos via detection of Cz^B. 

To illustrate the impact of the existence of relic right-handed neutrinos, we draw the expected spectrum 
of electrons by assuming that the background arises only from the beta decays of tritium. Moreover, to 
account for the effects due to a finite energy resolution, we follow the approach given in Ref. m by 
convolving both the beta-decay spectrum and the true CzzB signal with a Gaussian envelope of FWHM A. 
Under these assumptions, we consider two benchmark scenarios by taking (AjTtzq) as (0.15 eV, 0.25 eV) 
and (0.05 eV,0.07 eV), respectively. These two sets of benchmark values are chosen to respect the rule 
of A < 0.7mQ, which was found to be the necessary condition of discovering relic neutrinos from beta- 
decay background M- In addition, for the first benchmark scenario we choose its energy resolution to be 
compatible with the PTOLEMY experiment, while in the second one we reduce the energy resolution so 
as to meet the neutrino mass requirement S < 0.23 eV at the 95% C.L. from the latest Planck result [19] . 

The expected spectra of electrons for the above benchmark scenarios are drawn in Fig. [21 given different 
values of the total capture rate. Here, the differential capture rate dP/dEe has been shown as a function 
of the kinetic energy of electrons calibrated by which is the beta-decay endpoint in the limit of 

massless neutrinos. The standard scenario of massive Dirac neutrinos with T^ = 4 yr“^ is represented 
by the red solid curve, while the new-physics scenario with r§ = 5.1 yr“^ by a blue solid curve. For 
comparison, the capture rate Pjy; = 2rQ in the case of Majorana neutrinos has also been given as a green 
solid curve. As one can see, because of the enhancement of the capture rate, the observed spectrum in 
the latter case has a peak higher than the former one. Therefore, if one can experimentally resolve such 
a difference in the spectrum, a discrimination between the standard Dirac scenario and its modification 
might be possible. We next study such a discrimination quantitatively by taking the first benchmark 
scenario as an example, as it more closely resembles the proposed PTOLEMY experiment. First, we take 
the capture rate in the standard Dirac case (i.e., T = 4yr“^) as the true signal rate, and generate the 
true spectrum of electrons for a specific choice of data-taking period. Then, we fit this true spectrum by 
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Figure 2: Expected spectra of emitted electrons in PTOLEMY with two sets of energy resolution A and 
neutrino mass tuq. The black dashed curves are the background contribution from beta decays, while 
the colored dashed curves are the possible signals with three different capture rates assumed, i.e., red 
for the standard Dirac case T^ = 4 yr“^, blue for the case with the inclusion of right-handed neutrinos 
Tq = 5.1 yr“^, and green for that twice the standard Dirac contribution. The observed spectra are denoted 
by colored solid curves, which are the sum of the background and signal contributions. 


two parameters, namely, the capture rate T and neutrino mass uiq. In the fitting, we consider a region 
of interest that spans from the zero of and towards the signal end, and take as many bins 

as possible to sufficiently cover the signal region, where a bin size of 0.15 eV identical to the energy 
resolution has been chosen. Therefore, for the first benchmark scenario, the energy region of interest is 
[0, 0.45] eV, divided into three bins. It is worthwhile to point out that this method is independent of 
the location of signal peak (i.e., the neutrino mass). The number of events in each bin is then used to 
calculate the probability distribution of the fitted capture rate T and neutrino mass ttt-q, by assuming a 
Poisson distribution of the event number. 

By carefully inspecting the two-dimensional probability distribution of P and tuq, we find that the 
neutrino mass can be very precisely determined. In Fig. [Sj we plot the normalized probability distribution 
of the neutrino mass rug by marginalizing over the capture rate F. As one can see, even with one year 
of data-taking, the fitted rug is well peaked at the true value of rug = 0.25 eV and its width is extremely 
small. This is due to the fact that if the fitted neutrino mass slightly differs from the true value, the 
corresponding spectrum, especially the beta decay part, would have a drastic distortion, resulting in a 
good discriminating power for the neutrino mass. However, one should also keep in mind that in reality 
there might be large systematic uncertainties in modeling the beta-decay spectrum, such as the uncertainty 
of the energy resolution. 

Having reconstructed the neutrino mass from the spectrum of electrons, we now examine the discrim¬ 
ination among different capture rates. To do so, we first derive the probability distribution of the capture 
rate F by marginalizing over the neutrino mass rug. Then, for a given capture rate F, the p-value is 
calculated as the exclusion probability by integrating the distribution function over the values larger than 
F [48]. In Fig. [4| we show the p-value as a function of the capture rate. Several different data-taking 
periods are considered. It can be seen that for the non-standard Dirac scenario with F^ = 5.1 yr“^ (the 
vertical dashed line), running PTOLEMY-like experiments for about 5 years is only able to distinguish it 


II 















Figure 3: The probability distribution of the neutrino mass itiq, where the true value uiq = 0.25 eV and 
one year of data have been assumed and the capture rate F has been marginalized over. 



r(yr-') 


Figure 4: The p-value is shown for the test value of the capture rate F, where the true signal has been 
assumed to be Fd =4 yr“^ in the standard Dirac case. The energy resolution A and the true value of 
neutrino mass uiq are chosen to be 0.15 eV and 0.25 eV, respectively. The vertical dashed line is for the 
reference value F^ = 5.1 yr“^. 


from the standard one at la level, and to reach a 3a exclusion a detection time of even 20 years is not 
enough. Therefore, one may consider increasing the target mass in order to reduce the required data-taking 
time. As for the discrimination between the Dirac and Majorana scenarios, from Fig. HI we observe that a 
data-taking time of 5 years is sufficient for distinguishing these two scenarios at 3 ct level, while a detection 
time of about 10 years is required to achieve a 5a significance. 

Although our statistical analysis is performed for a particular choice of energy resolution and neutrino 
mass, it may apply to other similar scenarios, where the signal spectrum is well separated from the beta- 
decay background (i.e., there exist energy bins in which signal contributions are dominant). However, for 
the cases where the signal and background are not easy to separate, one would need a longer data-taking 
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time. In this sense, the results given in Fig. 0] may be regarded as the most optimistic case, and the 
indicated detection time can be viewed as lower bounds for more general choices of energy resolution and 
neutrino mass. 

In the above discussions, we have used the average number densities of neutrinos to calculate the 
capture rates of Ci^B, but it should be noticed that an overabundance of neutrinos within the dark matter 
halo is possible via the gravitational clustering [491114j . Depending on the absolute neutrino masses, the 
clustering effects could be significant. For the absolute scale of neutrino masses ttiq = 0.15 eV, the capture 
rate will be enhanced by a factor of 1.4 if the Navarro-Frenk-White profile [50] of the dark matter in our 
galaxy is assumed, or by a factor of 1.6 for the Milky Way model m- Hence, the unknown dark matter 
profile leads to the remarkable uncertainty of gravitational clustering, which can be even larger for heavier 
neutrinos, rendering the discrimination between Dirac and Majorana neutrinos extremely challenging. 
Additionally, in the presence of relic right-handed neutrinos, the difference between the capture rates Fj^ 
and F]y[ becomes even smaller. Therefore, the uncertainty in gravitational clustering effects and a possible 
cosmic background of right-handed neutrinos will diminish the experimental discriminating power on the 
nature of neutrinos. 

Finally, it is also interesting to think about a possible way to discriminate the thermal and non-thermal 
production mechanisms of right-handed Dirac neutrinos. The neutrino energy spectra are quite different 
in these two cases. For instance, the non-thermal distribution proposed in Ref. m favors low-energy 
neutrino states, implying the importance of measuring the velocities of hnal-state electrons. However, the 
terms involving the incident neutrino momenta in the capture rate are neglected, as they are typically very 
small and further suppressed by the heavy mass of the target nuclei. Hence, the final capture rate F^^^g 
mainly depends on the number density of the incoming relic neutrinos, while the velocity distribution or 
the detailed energy spectrum of relic neutrinos is almost irrelevant. 

6 Summary 

In this paper, assuming massive Dirac neutrinos, we have considered the impact of right-handed neutrinos 
on the detection of cosmic neutrino background in the future PTOLEMY experiment, in which the beta- 
decaying tritium will be used to capture background neutrinos. 

First, we demonstrate that the production rate of right-handed neutrinos is extremely small in the 
extended SM, as already found in Ref. |20|. Although our calculations show that the rate can actually be 
larger by several orders of magnitude than that in Ref. |20| , it is still insufficient for right-handed neutrinos 
to be thermally produced. Second, in assumption of thermal right-handed neutrinos present in the early 
Universe, we find that the cosmological constraints on the effective number of neutrinos require them 
to decouple from thermal bath at latest in the epoch of QCD phase transition. When the cosmological 
upper bound < 0.53 is saturated, we obtain the right-to-left ratio of neutrino number densities 

~ 0.28. Namely, the relic density of right-handed neutrinos z/g and left-handed antineutrinos 
Eg can be as large as 95 cm~^. Third, we present two possible scenarios to realize thermal production 
of right-handed Dirac neutrinos. The first one is just to utilize the magnetic dipole moments of massive 
Dirac neutrinos in the SM, given primordial magnetic fields generated in the electroweak phase transition. 
The second one is to introduce secret interactions among right-handed neutrinos and cold dark matter, 
which could help solve the small-scale structure problems. Finally, we examine how the presence of right- 
handed neutrinos affects the capture rate of cosmic neutrino background. Quantitatively, for massive Dirac 
neutrinos, the capture rate can be enhanced from Fj-, = 4 yr“^ to 5.1 yr“^ in the PTOLEMY experiment 
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with 100 grams of tritium [Mils]. To observe the impact of right-handed neutrinos, a data-taking time 
about 5 years is needed to reach a statistical significance of la in the most optimistic case. Therefore, it 
seems difficult to exclude or prove the existence of relic right-handed neutrinos in the near future. 

The proposal of PTOLEMY experiment and its great physics potential have stimulated us to take 
more seriously the detection of cosmic neutrino background and some related issues, such as the presence 
of right-handed Dirac neutrinos and the discrimination between Dirac and Majorana neutrinos. Further 
progress in both experimental and theoretical studies in this direction will hopefully extend greatly our 
knowledge on the intrinsic properties of neutrinos. 


Note Added: 

When the present paper was in preparation. Ref. m appeared in the preprint archive. Although both 
our paper and Ref. m have considered the impact of relic right-handed neutrinos on the detection of 
cosmic neutrino background, the production mechanisms for right-handed neutrinos are different and 
complementary. After finishing this work, we became aware of Refs. [571[5S1[22], where the right-handed 
neutrinos as dark radiation have been discussed in a U{1) gauge model with a heavy Z' gauge boson. 
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